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Abstract. Using the input energy per cycle as a quantifier of stochastic resonance (SR), we show that 
SR is observed in superharmonic (hard) potentials. However, it is not observed in subharmonic (soft) 
potentials, even though the potential is bistable. These results are consistent with recent observations 
based on amplitude of average position as a quantifier. In both soft and hard potentials, we observe 
resonance phenomenon as a function of the driving frequency. The nature of probability distributions of 
average work are qualitatively different for soft and hard potentials. 
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1 Introduction 

Stochastic resonance is an exclusively nonlinear phenomenon 
where the combined effect of the noise and the nonlin- 
earity gives rise to an enhanced response of the system 
at a particular frequency of an external periodic drive. 
It has been found to be of fundamental importance not 
only in physics PJ(I] but also in biological systems, from 
the mechanoreceptor cells in crayfish to the functioning of 
sensory neurons in humans. This is in sharp contrast to 
the general trend of a noise to cause the effect of a signal 
to fade. A typical model used to study this behaviour con- 
sists of a bistable potential in which a Brownian particle is 
present. The particle is in contact with a thermal bath of 
temperature T. This system is driven by a periodic drive 
with a given frequency, f(t) = Asinivt. Now, the initial 
system without the drive has an intrinsic escape rate of 
going from one minimum of the potential (V(x)) to the 
other. This is given by the Kramers escape rate [5J[S] : 
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where /3 = 1/ksT, fcg being the Boltzmann Constant, 
C is a constant that depends on the system parameters, 
and AV is the barrier height (height difference between 
the minimum and the maximum of V(xj). The escape 
time will then be given by the inverse of the escape rate: 
As the external periodic drive is switched on, 
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1 K 



in general its time period t w will not be in synchroniza- 
tion with the escape time of the particle over the barrier. 
However, if the noise strength or temperature is varied, at 
a certain value of temperature, t w will be exactly equal 
to 2tk- Now we will have proper synchronization of the 
dynamics: when the right well becomes deeper compared 
to the left well, the particle hops into the right well with 



a high probability, and vice versa. Under this condition, 
the system absorbs maximum energy from the drive. This 
phenomenon is termed as stochastic resonance (SR). Vari- 
ous quantifiers of SR have been proposed in literature [1] : 
signal-to- noise ratio (SNR), hysteresis loop area (HLA), 
spectral power amplification (SPA), position amplitude of 
the particle (aj), phase lag of the response with the ex- 
ternal drive (<fi), etc. In [7j, a relation has been derived 
between SNR, HLA and SPA which is given by: 



SNR x HLA = - 
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The above relation shows that some of the above quan- 
tifiers are related to each other. We observe that a simi- 
lar relation is present that connects the quantifiers mean 
thermodynamic work, x and <j), and as a result the three 
cannot act as independent quantifiers of SR. 

In the present work we use the mean input energy per 
drive period as a quantifier for SR [8HTTj. As a function 
of noise strength, we observe suppression of SR in soft po- 
tentials and distinct peak signifying SR in hard potentials. 

Since it has been established that SR is a bonafide 
resonance it should show maximum in the quantifiers as a 
function of drive frequency as well. Interestingly, whereas 
the mean input energy per period, (W), shows peaking 
behaviour, the average amplitude, x does not. As opposed 
to the behaviour of (W) as a function of noise strength, 
for soft and hard potentials, its behaviour with frequency 
is same for both kinds of potentials. In particular, (W) 
shows peaking behaviour and x decreases monotonically 
as a function of frequency for both hard and soft poten- 
tials. We further investigate the probability distributions 
of work for various temperatures jT0J|Tl] , for both hard 
and soft potentials, and point out the qualitative differ- 
ences between the nature of these distributions. 
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Fig. 1. The shapes of the potential for different values of q, 
with k = 0.2. 



2 The system 

We consider a Brownian particle described by the over- 
damped Langevin equation: 



1 x = -V'(x) + f(t)+m, 



(3) 



where V'(x) = dV g ^ , £(t) is Gaussian distributed white 
noise having the properties = and (£(£)£(*')) = 

2jkBTS(t — t'), 7 being the viscous drag in the medium. 
The strength of the noise is described by the thermal 
energy ksT. The expression for the bistable potential is 
given by [12] 

2 \r\ q 

V(x)=e~ x +k l -j-, (4) 

where the parameter k has been set equal to 0.2 through- 
out the manuscript, which sets the barrier height at ap- 
proximately AV = 0.67 (all variables are dimensionless). 
The parameter q is used to modulate the steepness of the 
walls of the potential. In other words, as q is increased, 
the slope of the potential wall increases as shown in figure 
[TJ so that the particle is more confined in-between the two 
minima. The external drive is given by f(t) = Asinwi. 

If we plot the amplitude of mean position as a function 
of temperature, i.e., x(T) vs T, then the resonance gets 
suppressed when q < 2. Analytically, this follows from the 
following approximate expression for (x) in the nonequi- 
librium steady state (where the intrawell dynamics has 
been ignored) [1]: 



where 
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Here, the angular brackets (• • •) represent ensemble av- 
eraging over a large number of phase space trajectories. 
(x 2 )q is the variance of the position of the particle in ab- 
sence of any drive, i.e., subjected to the unperturbed po- 
tential, m is the Kramers escape rate whose expression is 



given by 0E] 



where ±x rn and are the positions of the minima and of 
the maximum of the potential, AV is the barrier height, 
and V" is double derivative of the potential function with 
respect to x. 

It can then be shown from ^ that as T — >• oo, the 
behaviour of the amplitude of mean position is given by 



lim x(T) - T 



2/9-1 



(9) 



Now it can easily be seen that if q > 2, then x(T) goes 
to zero for large T. Of course, as T — > 0, x — > as well, 
because the particle hardly deviates from its equilibrium 
position (for details refer to [H]). This means that there 
must be a maximum in-between these two limits - a sig- 
nal for resonance. However, if q < 2, the particle travels 
large distances away from the minima, so that the x grows 
monotonically with temperature, and stochastic resonance 
is not observed. The case q = 2 is the marginal case. 

Following stochastic energetics [16] , the thermodynamic 
work done on the system is given by 



W 



T dV{x,t) 
dt 



dt 



X{ t)^>dt, 



(10) 



where V(x,t) = V(x) - xf(t). 

The average work done over time r is 

(H '> - - L' M, " d ^ 

Now, using the approximate expression for (x(t)) (eq. (O), 
we get 



(W) 
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Thus, we have arrived at a relation that connects three of 
the proposed quantifiers of SR: (W), x and (f>. Plugging 
in the expressions for x(T) (eq. ©) and 4>(T) (eq.©), we 
find 

27rA 2 (x 2 ) U r K 

Eq. (TT21 predicts SR as a function of temperature. 
Later on we use this to analyze our numerical results. 



3 Results and discussions 

The plots of average thermodynamic work done on the 
particle by the drive per period t u have been shown, for 
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A=0.1,v=0.02 





Fig. 2. Average work as a function of temperature for A = 0.1 
and v = 0.02 for different values of q. SR is observed for all 
values of q excepting q = 1.5. 



different values of the parameter q in figure [3] These plots 
have been obtained numerically by using the Heun's method 
[T3"] . We have ignored the initial transients and have eval- 
uated the work over many cycles (~ 10 5 ) using a single 
long trajectory of the particle. We find that the plots qual- 
itatively show the same features as shown by the position 
amplitude with temperature [12] . The parameters used 
have been given in the figure captions. At very low tem- 
perature, the particle can see only a single well, the barrier 
height being too large for it to cross. Thus, only intrawell 
dynamics is dominant under this condition and as a result 
the work done on the particle is very small and goes to 
zero as T — »■ 0. At the other extreme, T — > oo, the barrier 
becomes negligible compared to the thermal energy of the 
particle. Now the random motion of the particle becomes 
so large that the synchronization gets washed away. This 
happens only for strong confining potential with q > 2. 
Thus for the hard potentials, we get a clear resonance 
peak. Moreover, our numerical result shows that the tem- 
perature at which SR peak occurs is consistent with the 
condition r w = 2tr-. 

However, for values of q < 2, the particle travels far 
from the left (right) of left (right) minimum and takes a 
long time to return. So it is expected that the distribution 
of passage time above the barrier will be very broad, so 
that the mean passage time ceases to be a good variable, 
being dominated by a large dispersion. As a result, the 
synchronization condition of escape rate with drive pe- 
riod is never satisfied. Thus these plots do not show the 
characteristic maxima of SR. 

In figure [31 the numerically obtained plot for (W) vs T 
has been compared for the superharmonic potential with 
q = 6 with the following two analytical expressions for 
average work in steady state that are commonly used in 
the literature [T2"] : 



Fig. 3. Comparison of temperature dependence of (W) ob- 
tained numerically with the analytical results. (Wi) (the 
smooth solid line) is the expression obtained by taking into ac- 
count the temperature dependence of the variance in position 
of the particle, (x 2 )o- (W2) (the smooth dashed line) is the ex- 
pression for average work obtained with the variance replaced 
by ij,. The curve labelled (W) is the numerically generated 
one. 



In the expression for (Wa), we have replaced (x 2 )o by 
x 2 m , where x m is the position of the minimum of the poten- 
tial. However, from the figure we find that (Wi) matches 
the simulated curve reasonably well, whereas (W2) devi- 
ates by a larger extent. This can be understood as fol- 
lows [12] ■ In the derivation of the expression for (Wi), one 
assumes that we are dealing with a strictly two-state sys- 
tem, where the particle's position distribution is the sum- 
mation of two delta functions, x^ then becomes the vari- 
ance of the distribution having two 8- functions equidistant 
from the origin: (x m — 0) 2 = x 2 n . Evidently, at any finite 
temperature the above distribution will be incorrect, ow- 
ing to the softness of the double well potential. Thus, we 
need to incorporate into our expression the temperature 
dependence of the variance in particle position, which has 
been done in deriving the expression for (W\). 




(Wi) 



2ttA 2 (x 2 ) ur K 



T 
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(13a) 



Fig. 4. Phase difference <j> between drive and response, as a 
function of bath temperature. The peak in the curve becomes 
less distinct as the value of q is lowered. 
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Previously the phase lag <j> of the response with drive 
(eq. (0) has been used to detect stochastic resonance |15j . 
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The variation of <p with temperature has been shown in 
figure [4] We find that the curves for q = 2, 4 and 6 show 
prominent maxima, whereas the curve for q = 1.5 mono- 
tonically increases from zero and then saturates to an up- 
per limit. Systems exhibiting SR show a peak in the phase 
lag cf> as a function of noise strength. However, the opti- 
mum value of noise intensity at which peak occurs does 
not coincide with SR peak for different quantifiers. The 
bell-shaped dependence reflects the competition between 
hopping and intrawell dynamics pQ. However, for q = 1.5, 
we do not see a bell-shaped curve, implying no clear-cut 
time scale separation between hopping and intrawell mo- 
tion. 



3=2 



Fig. 5. Average work as a function of frequency for A = 0.1 
and T = 0.3 for different values of q. SR is observed for all 
values of q. 



n 




Fig. 6. position amplitude as a function of frequency for dif- 
ferent values of q. The monotonic decrease in x is apparent. 



We now study the SR quantifiers (W) and func- 
tion of frequency of drive. The major difference between 
x and (W) as quantifiers is observed as a function of fre- 
quency of external drive. In figure [5] we have plotted the 
mean work versus driving frequency. Whereas (W) shows 
a peak for all values of q (figure [5]), x decreases monoton- 
ically with frequency for any value of q (see eq. ([6])), as 
can be seen in figure [51 It may be emphasized here that 
the general trends of (W) and i as a function of w do not 



depend on the nature of the confining potential. This is 
in contrast to the behaviour of (W) and function 
of temperature, which crucially depends on the softness of 
V(x). 

Computing the derivative of (W) with respect to uj 
from the expression (|12[) . we find the maximum to occur 
precisely at uj = 2tk'- 



d(W) 



dui 







LO r , 



2r K - 



(14) 



However, the proper condition for synchronization is 

t w = 2t k => u>sb. = nr K . (15) 

On calculating the escape rate for the superharmonic (q = 
4) potential at T = 0.3, we find tk ~ 0.03. Figure [5] shows 
the peak to occur at u) ~ 0.06 ks 2r^, consistent with eq. 




Fig. 7. Comparison of frequency dependence of (W) obtained 
numerically with the analytical results. (Wi) (the smooth solid 
line) is the expression obtained by taking into account the vari- 
ance in position of the particle, (a; 2 )o, at the fixed temperature 
T = 0.3. (Wz) (the smooth dashed line) is the expression for 
average work obtained with the variance replaced by x^. The 
curve labelled (W) is the numerically generated one. 

In figure [71 the numerically obtained plot for (W) vs 
u) is compared with the analytical expressions (|13al) and 
(|13bl) . For this we have used the superharmonic potential 
with q = 6. We observe that the analytical expression 
using eq. (|13ap fits better compared to (|13b[) . the reason 
being the same as explained earlier in the context of figure 
Now we turn our attention to the power applied to 
the system. 

In figure [HI we have plotted the power applied to the 
system by the drive versus the drive frequency. We find 
that for all values of q, the curves are monotonically in- 
creasing. From figure [2] (scaling the y-axis by the constant 
parameter w), we observe that the average power exhibits 
peak for q > 2 as a function of temperature. Thus, power 
cannot be used as a quantifier of bona fide SR [T3] for 
q > 2. 

In figure [9l we have plotted the real time trajectory of 
the particle at T — 0.3 (around SR for q=4). In (a), be- 
cause of subharmonic potential, the particle travels large 
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Fig. 8. Average power as a function of frequency for A = 
0.1 and T = 0.3 for different values of q. The curves are all 
monotonically increasing with frequency. 




(a) 




(b) 

Fig. 9. (a) Position as a function of time for a given trajectory 
is shown, for the q — 1.5 potential. We find that the particle 
travels large distances from either minima due to softer con- 
finement, (b) Similar plot for q = 4. The hard confinement 
effectively contains the particle in a smaller region, and the 
motion is more synchronized. 
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Fig. 10. Work distributions at different temperatures for q=4. 
Other parameters are A = 0.1 and v = 0.02. 



distances away from the minima and spends more time in 
the wings of the potential (x > x m or x < —x m ) over a du- 
ration of many cycles of the applied force without passing 
over the barrier. This is clear from the figure. Thus, the 
question of synchronization between the applied force and 
particle hopping does not arise, hence the absence of SR. 
On the other hand, in (b), the superharmonic potential 



(q = 4) helps in more efficient confinement of the particle 
so that the proper synchronization between the drive and 
the particle trajectory is attained. 

In figure HOI we show how the work distribution P(W) 
changes as a function of T for the superharmonic potential 
with q = 4. At small temperature (T = 0.1, figure [TP] (a)), 
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playing dominant role and hence the distribution becomes 
broader (figures [10] (b) and (c)). Additional peak appears 
towards right mainly due to the interwell motion. A third 
peak also appears in the negative side. For large values of 
temperature beyond SR point, the dynamics is dominated 
by interwell motion and P(W) tends towards a Gaussian 
distribution (figure [10] (d)). 

The probability distribution for work has finite weight 
for negative values of W . These negative values correspond 
to the trajectories where the particle moves against the 
perturbing ac field over a cycle. The existence of finite 
weight for negative work values is essential to satisfy the 
recently discovered fluctuation theorems. 

In figure [TT1 we have plotted the variations in the work 
distribution with temperature, for the subharmonic po- 
tential (q = 1.5). This time, however, we do not find the 
appearance of prominent double and multiple peaked dis- 
tribution as was observed for the superharmonic potential. 
This is because the particles travel higher distances during 
intrawell as well as interwell motion, which gives rise to 
work values varying over a wide range. Moreover there is 
no clear-cut time scale separation between intrawell and 
interwell motion. Thus the two distinct peaks that were 
observed for q = 4 have got merged in q = 1.5 case. 
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4 Discussion and conclusions 

In this paper, we have studied the phenomenon of stochas- 
tic resonance in a bistable potential, using the mean input 
energy per cycle (or the mean work done per cycle) (W) 
as a quantifier of resonance. We find that the system ex- 
hibits SR as a function of temperature for q > 2, but 
does not show SR for subharmonic potentials. This be- 
haviour is further verified by studying the phase lag <f>. 
Thus bistability is necessary but not sufficient condition 
for the observation of stochastic resonance. This result is 
consistent with the findings in 12 . However, in both the 
superharmonic and subharmonic potentials, the work ex- 
hibits resonance peak whereas the average amplitude of 
mean position decreases monotonically as a function of 
frequency. This is quite different from the trends of (W) 
and X ciS el function of temperature, which is sensitive to 
the nature of the confining potential. We have shown that 
the average power delivered to the system is not a good 
quantifier for bonafide resonance [14) . Our further investi- 
gation reveals qualitative differences in the nature of dis- 
tributions for hard and soft potentials. 



Fig. 11. Work distributions at different temperatures for 
q=1.5. Other parameters are A = 0.1 and v = 0.02. 



the particle sees only a single well, the barrier height being 
too large for it to cross. Thus the work done is entirely 
due to intrawell dynamics and the distribution is almost 
Gaussian. Occasional excursion of the particle into the 
other well is clearly reflected as a small hump at higher 
values of W . As T increases, interwell dynamics starts 
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